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A Unified Theory of Dynamical Mean Field and Spin Fluctuations 

Tetsuro SasoQ 

Department of Physics, Faculty of Sciences, Saitama University, Urawa 338-8570 
(Received May 25, 1999 ) 

A new approach is proposed which encompasses the dynamical mean field theory (DMFT) 
for strongly correlated electron systems and the self-consistent renormalization (SCR) theory of 
spin fluctuations. The latter is incorporated into DMFT as a 0(l/d) correction, where d is the 
spatial dimensionality, but in a phenomenological way. The DMFT is completely recovered when 
the effect of spin fluctuation is omitted, whereas the properties at the quantum critical point 
(QCP) is the same as the SCR theory. An example of calculation is presented for the single 
band Hubbard model and the non-Fermi liquid behavior is demonstrated at QCP. 
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§1. Introduction 

It is one of the most important tasks in condensed 
matter physics to construct a theory which can describe 
the electronic states and various physical properties of 
the strongly correlated electron systems (SCES) includ- 
ing high temperature superconductors, heavy electron 
systems, organic conductors and so on. Various numer- 
ical methods so far used, e.g., numerical renormaliza- 
tion group, density matrix renormalization and quantum 
Monte Carlo methods, are valid only up to 0, 1 and 2 
spatial dimensions, respectively, and extension to three 
dimensions or to realistic systems with complicated band 
structure seems difficult. _ 

The dynamical mean field theoryEJ 1 (DMFT) is a most 
powerful method among such theories for describing the 
effect of strong correlation. It becomes exact in the limit 
of large spatial dimension d — > oo, where the lattice prob- 
lem is reduced to solving an impurity problem embedded 
in an effective medium self-consistently. It can be re- 
garded as a natural extension of the mean field theory 
for spin systems to the itinerant electron systems and 
offers a theoretical description of the latter beyond the 
Hartree-Fock approximation. Most properties of Ihrcc 
dimensional systems, including Kondo insulatarspEP are 
expected to be described well by this method.El 

DMFT, however, neglects corrections due to finite d, 
which are the fluctuations from the mean field and may 
become important in the vicinity of the phase transition 
points. Several attempts to extend. DMFT have been 
made so far. Schiller and Ingersentffi included 1/d cor- 
rection perturbatively. Hettler, et ala considered a clus- 
ter embedded in an effective medium. These approaches, 
however, do not seem to have reached a level of practical 
use. 

In the case of the spin fluctuation, the_so-called self- 
consistent renormalization (SCR) theoryEP is proposed, 
first for d-electron systems, where the correlation is 
not very strong. SCR theory can be microscopically 
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derived, 13 Ef but the recent phenomenological construc- 
tion© described with only small number of physical 
parameters, seems much simpler and having wider ap- 
plicability. It was later extended to the f-electron sys- 
tems where the correlation is much stronger .Et This was 
done also on a phenomenological basis. Recently, much 
attention^ has been paid to the quantum critical point 
(QCP),E£P at which, e.g., antiferromagnetic long range 
order is destroyed by doping, applying pressure, etc. It 
has been proved that SCR can describe physics around 
QCP successfully.tlr On the other hand, since SCR de- 
scribes a system by only small number of parameters, 
details of the structure specific to the system may not 
be taken into account. Miyake and NarikiycO extended 
SCR to include the effect of nesting and the sharp struc- 
ture of the density of states. It would be valuable if SCR 
theory for SCES could be rebuilt on a more microscopic 
basis. Purpose of the present paper is to propose such 
a theory which starts from DMFT and is extended to 
include SCR, or that which interpolates between DMFT 
and SCR. _ 

Kuramoto and MiyakaiiJ' proposed a phenomenologi- 
cal theory of heavy electron systems which emphasizes 
the dual character of the f-electrons and separates the 
mixture of conduction and f electrons into the localized 
moment part and the heavy quasi-particle part, which 
are interacting with each other. Extending the Fermi liq- 
uid theory, they constructed the dynamical susceptibil- 
ity of the impurity Anderson model, and then proceeded 
to the duality description of SCES. Basic assumption 
of Moriya and TakimotoS in extending SCR to SCES 
stands on ±he_above theory. 

OhkawaoliS* discussed an extension of DMFT to in- 
clude intersite spin fluctuations through 0(1 /d) correc- 
tion in detail. His formulation, however, is rather compli- 
cated and needs further approximation or simplification 
in actual calculations. 

Therefore, we propose in this paper a simpler scheme 
which unifies DMFT and SCR and interpolates between 
them. DMFT-js treated by the iterative perturbation 
theory (IPT),!^ which is simple but valid for wide pa- 
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rameter range of SCES. The phenomenological construc- 
tion of SCR, to which its simplicity owes, is kept as much 
as possible in the present formulation. 

The present paper is organized as follows. In §2, the 
methods of DMFT and IPT are summarized. In §3, a 
modified IPT scheme, which can be used for electron-hole 
asymmetric case, is explained. In §4, a unified theory of 
DMFT and SCR theory is presented. Numerical calcu- 
lation for the Hubbard model is presented as an example 
of the calculation. Finally, discussion and conclusion are 
given in §5. 

§2. Iterative Perturbation Theory 

In the present paper, we investigate the single-band 
Hubbard model, 

ia i 

(2.1) 



no- 



where Sriirj — m a — {iiia) and E d = E d + U{n- a ) de- 
notes the Hartree-Fock level of d electrons. m a denotes 
the electron number at site i with spin a and (• • •) the av- 
erage. We do not apply magnetic field, so that {n irT ) does 
not depend on a nor i. (• • •) will be omitted. In DMFT, 
the hopping integral between the nearest neighbor sites 
is usually scaled as Uj — t/2vd. If we use the simple 
hyper-cubic lattice, it is well known that the density of 
states is reduced to the Gaussian form of the width t 
in the limit d — > oo. Furthermore, the lattice problem 
is reduced to solving the impurity problem embedded in 
an effective medium self-consistently. Here, the impu- 
rity can be regarded as being placed at the cavity site 
in the center of the effective medium. Its unperturbed 
state is described by the cavity Green's function G(e), 
the Green's function at a cavity site (interaction U is 
artificially removed there), which is related to the local 
Green's function G(e) as 



G(e) = (G( e )- 1 + S(e))- 



and G(e) is defined by 



G( £ )4r — — • 



(2.2) 



(2.3) 



N is the number of lattice sites and e k denotes the 
energy of conduction electrons. E(e) denotes the self- 
energy which describes the effective medium and is self- 
consistently generated at the cavity site by U. 
In IPT, the self-energy S(e) is calculated up to the 




Fig. 1. The second-order self-energy. 



second order of U (Fig.l) as a function of the Matsubara 



frequency as 

£( 2 >(ie) = C/ 2 T^G(ie + icj)n(kj), (2.4) 

UJ 

where e = (2n + 1)tvT, uj = 2mtT (n ^integer) and 

fi(iu;) = ~T^2 G(ie')G(ie' + iw). (2.5) 

e' 

The self-energy is analytically continuated to the real 
frequency as 

S (2) (e) = U 2 J dei J de 2 J de 3 p( ei )p(6 2 )p(e 3 ) 

x /(-ei)/(-e 2 )/(63) + /(ei)/(e 2 )/(-e 3 ) 
e - ei - e 2 + e 3 + i6 

where p(e) = (— l/7r)ImG(e + id) is the density of states 
of the cavity, 5 a small positive number and /(e) denotes 
the Fermi function. 

This equation is converted into the following form: 

/>oo 

£ (2) (e) = -kV 2 / dt e iet [(3(t)a{-t)P(t) 
Jo 

+ a(t)/3(-t)a(t)\ , (2.7) 
by introducing the transformation of /5(e): 

til } = /l d£ e " i£t/5(e)/(±e) - (2 - 8) 

In this paper, we use the semi-circular density of states 
p(e) = {2/ttW)^J\ — (e/W) 2 for the unperturbed elec- 
trons. Then the Green's function, eq.(2.2), is calculated 

as 

z + WW 2 - z 2 
E d — S(e). Combining with eq.(^.2]) we 



with z = e 
obtain 



G(e) = 



^G(e)' 



(2.10) 



These equations should be calculated self-consistently. 

§3. Modified IPT 

IPT reproduces the DMFT calculation via exact di- 
agonalization of the effective impurity in the case of 
electron-hole symmetry rather well,EZP in spite of the use 
of the second order perturbation. A reason is firstly that 
the second-order perturbation for the impurity Ander- 
son model reproduces the atomic limit correctly in the 
symmetric caseo-* Thus it succeeded in describing the 
Mott transition- of the Hubbard model on Bethe lattice 
at U /W w 3.113' Recently, it was pointed out that an 
another reason may be that IPT becomes exact _a[hen 
the width of the conduction band becomes zeroO In 
the electron-hole asymmetric cases, however, the second- 
order perturbation does not reproduce the atomic limit, 
and the positions of the Hubbard bands are wrong. 

To remedy .this point, a phenomenological interpola- 
tion schcmcL3Eil'E3'E3 was proposed to encompass the 
perturbative regime and the atomic limit. In this treat- 
ment, the following form is assumed for the modified 
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self-energy: 



E(e) 



l-BS( 2 )(e) : 



(3.1) 



where the parameters A and B are determined so as to 
reproduce the correct |e| — > oo limit and the atomic limit, 
respectively, yielding 



.4 



B 



n_ CT (l - n_ CT ) 



1 - 2n_„ 



Un-a(l 



(3.2) 



(3.3) 



Here, h a is the electron number calculated from p(e) in- 



stead of p(e) = — (l/7r)ImG(e 
function G(k, e) = (e — E d — 
of the atomic limit 

Gate) 



W). Then, the Green's 
E(e)) -1 reduces to that 

(3.4) 



e-E d -U' 



e-E d 

when [/ > W, but does not satisfies the Luttinger sum 
rule because of the phenomenological modification of the 
self-energy. This point can be remedied by tuning E d 
as the. effective level so as to satisfy the Luttinger sum 
rule.EBEif We, however, adopt a more simplified scheme, 
in which the sum rule is satisfied only approximately. 
Nevertheless, physically reasonable behavior is obtained 
for a wide range of parameters without tedious deter- 
mination of the effect hp^E^- This is achieved by sim- 
ply subtracting ReS(0)H ) at T = from S(e) at each 
temperature. We show the density of states for various 
values of E d with W = I and {7 = 2 in Fig. 2. The 
Fermi energy is set equal to 0. It is seen that almost 
correct shapes of the spectra are obtained for all cases. 
The Luttinger sum rule is stated as = n a , if we define 
the electron number by 



= ^0(/i-£ d -E(M)-efc), 




Fig. 2. The densities of states calculated by the modified IPT 
scheme for W = 1, U = 2 and various values of E d are displayed. 
The Fermi energy is set equal to 0. 



(3.5) 



where [i denotes the chemical potential. In Fig. 3, ?v in 
the modified IPT and in the Hartree-Fock approximation 
(HFA) are plotted as functions of E d . Since S(/x) = in 
our treatment, is approximately equal to n a in HFA, 
whereas n a ~ n a (HFA) in Fig. 3. Therefore, the Lut- 
tinger sum rule is approximately fulfilled in the present 
treatment. 

§4. Inclusion of Spin Fluctuations 

In this section, we first discuss the dynamical suscep- 
tibility of SCES. In extending ,SCR theory to the case 
of SCES, Moriya and TakimotoEf assumed the following 
form for the dynamical susceptibility, motivated by the 
duality picture of Kuramoto and MiyakeJij-' 



X (Q + q,to) = [xlM" 1 - J Q (T) + Aq 2 



(4.1) 



where xl( w ) denotes the local susceptibility of each mag- 
netic ion and Q an ordering vector. In SCR, xl( w ) 1 i s 
expanded up to the linear term of u>, yielding 

XL 



x! CR M = 



i-i^AV 



(4.2) 




Fig. 3. Number of electrons as a function of E d calculated by the 
Hartrcc-Fock approximation and by the modified IPT scheme at 
T = 0. 
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below certain cutoff frequency u) c . Here, Tl is of the 
order of the Kondo temperature of the magnetic ion. 

Kuramoto and MiyaketS derived an approximate form 
of the dynamical susceptibility for the impurity Ander- 
son model on the basis of the Fermi liquid theory. We 
calculate it here in a similar but slightly different manner 
as 

n( w ) 



i -t(uj)u(uj)' 



(4.3) 



where II (w) is defined in the same way as II(a;) in eq.( |2.5| ) 
but with G(e) instead of G(ui). The Feynman diagram is 
shown in Fig. 4. T(u>) is the renormalized vertex function, 
for which we adopt an approximate form, 

r < w > = irmry (4 - 4) 

which corresponds to taking fan type vertex corrections 
into account. This form of T(u>) has a useful property. 
Since T(lu) — > FI(u;) — 1 in U — > oo, the denominator of 
Xl(^) in cq.([Q|) does not vanish as far as U < oo and 
the magnetic instability is avoided. This is a necessary 
property when one wants to calculate x( w ) °f the single 
Anderson impurity for large U. Eq.( |4.3| ) can be rewritten 
as 



XL(w)=n(a;)[l + tai(w)], 



(4.5) 



which reproduces precisely the perturbation series ex- 
pansion of x(w = 0) up to the first order of U for the sym- 
metric and asymmetric cases of the Anderson model. E$ 




Fig. 4. Feynman diagram for XL ^ s snow where the bubble de- 
notes II and the square denotes T. 



In contrast to SCR theory, the present theory retains 
full dynamical structure of Xl( w )- But because of the 
approximate vertex correction, xl( w ) does not satisfy 
the sum rule, 



r-2 _ 3 
L - ~ 

7T 



dtj[l + 2n(u;)]ImxLM - S{S + 1), (4.6) 



where S — 1/2. In order to fulfill this sum rule, we 
modify xl( w ) as follows, 



XlM = [XL(w) 1 - iCuj]~ 



(4.7) 



with the parameter C(T) to be adjusted at each tem- 
perature. This may correspond to modifying Tl to an 
effective value to satisfy the sum rule. Here and hence- 
forth, we omit the prime on Xh- 

We show numerical results for xl( w ) at T = and 
Xl( w = 0, T) in Figs. 5 and 6, respectively. We have used 
the IPT calculation for — —0.5 in the preceeding 
section. For comparison, we also show H(u> — 0, T) in 



Fig. 6. It is seen that Xh{T) is enhanced over H(T) due 
to the vertex correction. 




Fig. 5. The real (solid line) and imaginary (broken line) parts of 



the local susceptibility xj j (w) are plotted for U 
and T = 0. 



2, E A 



-0.5 




Fig. 6. Temperature dependence of the local susceptibility 
XL (a; = 0,T) is plotted. TI(uj = 0,T) is also shown by the 
broken line for comparison. 



In DMFT, the dynamical susceptibility xil^) of 
the lattice can be obtained once the effective impurity 
problem was solved. lP The procedure, however, needs 
rather tedious calculations. OhkawaQ and Miyake and 
NarikiycO discussed a general form of x(g,w) and de- 
rived the following equation 



(4.8) 



where J(q,u>) — U 2 Axo(q, w) + A(q, ui) in the strong cor- 
relation limit. Axo(q, w) denotes the intersite part of the 
dynamical susceptibility of the lattice, 



Xo(<?,iw) 



-T^G(k,ie')G(k + q,ie' + hv), (4.9) 



Axo(?,iw) = Xo(«, V) 



N ^ 

q 



Xv{qM, (4.10) 



and X(q,ui) contribution from the mode-mode coupling 
of spin fluctuations. J(q,ui) includes both the mean 
field contribution to the exchange interaction, which is 
O(l/d ), and the fluctuations from the mean field, which 
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is 0(l/d). Here we assume that J(q, ui) can be expanded 
around a certain ordering vector Q as J(Q + q, u>) — 
Jq (T) — Aq 2 + ■ ■ ■ and w-dependence can be neglected. 
Then we obtain x(Q + 9S m the same form as eq.(l.l). 
We can determine Jq (T) by the condition that the spin- 
fluctuation amplitude 



Si 



dw[l + 2n(o))]IniXQ(w) 



(4.11) 



stays constant: (TV) = S*^(T), where Tjv is the or- 
dering temperature. Both quantum and thermal fluc- 
tuations are included in S*£(T). When there is no long 
range order, we set Tv = in the above condition. In 



eq.( lll ), Xq( w ) is defined by 



_3_ 



1 - yyoM tan 1 



(4.12) 



char- 



where y Q (u) = [xlM" 1 - Jq{T))/T a , T a = Aq% 
acterizes the stiffness of the spin fluctuation and qs is 
the wave vector at the Brillouin zone boundary. In the 
actual calculation we take q^ — 6tt 2 N/V, where V is the 
system volume. 

In Fig. 7, we show the temperature dependence of the 
parameter Jq(T) when Jq = Jq(T = 0) = 0.5 and 
0.7827= J c , which is a critical value and is equal to 
Xl (uj = 0,_T = 0)" 1 . J Q (T > 0) is determined by the 
condition (0) = 5£ (T) . 
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Fig. 8. The real (solid line) and imaginary (broken line) parts of 
the dynamical susceptibility XQ (<•<-') at T = are plotted. 
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Fig. 7. The temperature dependence of Jq (T > 0) is plotted for 
the cases Jq(T = 0) = 0.5 and J c . 



Fig. 9. The susceptibility xq(0) i s plotted as a function of T for 
Jq =0.5, 0.6, 0.7 and J c from upper to lower curves. 




Fig. 10. Feynman diagram for Xij ls shown, where the hatched 
bubble denotes xl an d the white bubble AXy ■ The filled circles 
denote U. 



The frequency dependence of Xq( w ) an( i the tempera- 
ture dependence of x(Q> 0) are plotted in Figs. 8 and 9 for 
Jq =0.5, 0.6, 0.7 and J c . It is seen that the low energy 
structure of xq(w) becomes singular and much narrower 
than xl( w )i an d xiQ-ify becomes proportional to T 3 / 2 
at QCP, as expected. u 

Once these calculations are done, we want to recalcu- 
late the self-energy. It must be done in a manner consis- 
tent with the IPT calculation. The dynamical suscepti- 
bility Xij ( w ) m the site representation may be expanded 



as 

xa = xtSij + u 2 xl^x%xl 

+ U^xl^xlxl^xl.xl + ■■ C/-\(4.13) 

k 

(Fig. 10), where x% denotes the irreducible susceptibility, 
X L = 1/U + xl an d the contributions from the succes- 
sive identical site are collected into xl = X%/{^ ~ Ux%), 
whereas the successive identical sites in the intermediate 
paths are excluded in the summation as indicated by the 
prime since &Xu = 0- The above equation can be solved 
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after Fourier transformation into 

xf» 



l-[/ 2 Ax (g,w)x L M' 



(4.14) 



, if 

XL 



which can be rewritten in the same form as eq.(4 
we define r J(g, ui) = [/ 2 Axo(<Z,w) and note that x£ 
in SCESM 

It may seem that one can calculate the self-energy of 
electrons as scattered by this spin fluctuation, but the 
calculation of the self-energy at the cavity site needs fur- 
ther modification. Namely, we divide eq.(4.8) into the 
form 



x(<2s w ) = xlH + 



1 - J(q,u)xh(^)'' 



(4.15) 



and replace xl(<-°) m the first term and in the numera- 
tor of the second term with n(w). Strictly speaking, the 
paths which go through the cavity site in the interme- 
diate steps must be subtracted, EZP but we neglect that 
effect in the present study. Then we obtain 



E(ie) 



U 2 T^G{ie 



(4.16) 



1 - 



ReE (J Q =J C ) 
ImE (J Q =J C ) 
Im£ (IPT) 








Fig. 11. The real (solid line) and imaginary (broken line) parts 
of the self-energy at T = is shown for Jq = J c . The imaginary 
part of the self-energy in IPT is also shown by the dash-dotted 
line for comparison. 



1 



q v ^ 



J{Q + q, iuj)fl 2 {iuj) 



■ i^)XL(i^) 



ri(iw) + [xq(iw) - XL(iw)] 



n(iw) 

XL(iw) 



(4.17) 



where we have replaced J{q,uj) with Jq(T) — Aq 2 . Note 
that the second term has the same criticality as xq(^)- 
If we set J(Q + q,ici>) = in this form, the self-energy 
eq.( 4.16 ) recovers the IPT result. 

For numerical calculation, the following expression is 
convenient (see Appendix): 



E(e) 
where 



lit) = - 

7T 



dt 



duj 



^(i)7(t)+o(i) 7 H)], (4.18) 



*n(o;)ImxQ(w + W). 



(4.19) 



The low energy regions of real and imaginary parts of 
S(e) are shown in Fig.ll, and the imaginary part is com- 
pared with that by IPT. It is seen that ImS(e) at e ~ 
becomes singular in contrast to the Fermi liquid behav- 
ior ImE(e) oc e 2 in IPT. Theoretical analysis leads to 
ImE(e) oc e 3 / 2 at QCP and T = 0. 

Using the above self-energy we calculate the density 
of states again, which is shown in Fig. 12. It is seen that 
the peak at Ep becomes thinner and singular due to the 
strong renormalization by the spin fluctuation. 

Finally, the total energy is calculated fromEJ 



1 



1 



E = 2 J de/(e) I -- ) Im[{e - -{Un- a + S(e))}G(e)]. 

(4.20) 

The specific heat is obtained by numerical derivative of 
£asC = dE/dT, which is plotted in Fdg.13. We see 
that C(T)/T cx -y/T at QCP as in SCR.H* 
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j Jq-J c 
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Fig. 12. The density of states at T = after inclusion of the spin 
fluctuation. The broken line is for Jq = 0.5 and solid line for 
Jq = Jc 




Fig. 13. Temperature dependence of C(T)/T is plotted for 
Jq =0.5, 0.6, 0.7 and J c . 
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The electrical resistivity, in the dimensionless form is 
calculated by the formulap 

/>oo 

R(T)= / dwn(w)[l + n(w)]Imxg(w). (4.21) 
Jo 

Since this formula is an approximate one, we do not take 
care whether we should use Xg( w ) in stead of xq( w )- 
Furthermore, the behaviors of xq( w ) an d Xq( w ) are the 
same at the critical point. The numerical result is shown 
in Fig. 14, which seems consistent with the theoretical 
analysis R(X) oc T 3 / 2 within numerical accuracy of the 
calculations' 



0.01 



H 

2 



0.005 




0.002 0.004 



0.006 



Fig. 14. Temperature dependence of the resistivity is plotted for 
Jq =0.5, 0.6, 0.7 and J c . 



§5. Discussion and Conclusion 

In this paper, we have proposed a new scheme which 
incorporates the dynamical mean field theory and the 
self-consistent renormalization theory of spin fluctua- 
tions. Starting from the IPT calculation (eqs. (2.7-10) 
and (3.1-3)) for DMFT, we calculated the renormal- 
ized local dynamical susceptibility xl( w ) approximately 
(eqs. (4.5-7)), and then constructed the dynamical sus- 
ceptibility xLojjjj) of the whole system in the .same form 
as Ohkawali-jtlif and Moriya and TakimotoB* Consid- 
ering the contribution of the mode-mode coupling in a 
phenomenological manner, we determined its tempera- 
ture dependence by the sum rule in parallel with the 
SCR theory (eqs. (4.11) and (4.12)). Both quantum and 
thermal spin fluctuations are taken into account. Since 
the theory is based on IPT, it is applicable to the strongly 
correlated systems. The self-energy of electrons due to 
the spin fluctuations is calculated by using the above 
x(q, w) with a correction at the cavity site in the mean 
field (eqs. (4. 17-19)), which ensures a consistency with 
IPT. The present theory encompasses the IPT calcula- 
tion for DMFT and the SCR theory. If we omit J{q,uj), 
we recover the IPT results, whereas the critical behav- 
ior at QCP is completely the same as SCR. Therefore, 
structures in the single-particle spectra specific to each 
material can be taken into account through the IPT cal- 
culation. Application to the Kondo insulators is also 
possible. Combination with band calculations will be an 
important subject in the future studies. When the nest- 
ing of the Fermi surface is strong, one has to take account 



of the frequency dependence of J(q, uj 

The present treatment still includes phenomenologi- 
cal constructions in incorporating the 1/d correction to 
DMFT. Furthermore, the self-energy is not calculated 
self-consistently, together with the effect of spin fluctua- 
tions. Namely, the 1/d correction was taken into account 
only in a perturbative manner. Recently, several at- 
tempts are done to systematically take account of the 1 jd 
corrections to DMFTrl'ElEZP Among them, Smith and 
SiEZP proposed an extended mean field theory which in- 
cludes the spin-spin coupling in a non-trivial way in addi- 
tion to the effective impurity action. Another important 
development in the effort to investigate the effect of fluc- 
tuatio^may be the two-particle self-consistent (TPSC) 
theoryj^r E£P which leads to a similar theory to SCR but 
can take account of both the charge and spin fluctua- 
tions. Combination of the present formulation with such 
theories may open a way to an extension of DMFT in a 
more microscopic fashion. 
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Appendix: Derivation of eq.(4.18) 



By analytic continuation, eq.( 4.16| ) is written as 



E(e) = U" 



duj , 



— [Pn( W )G R (e + W )ImX§H 



-f(e + W )ImG R (e + (A-l) 



r.A/ 



Using the spectral representation for G(e) and Xq( w ) 
we obtain 



E(e) = U 2 J duj J 



, n(u))+f(u) 
dv -p{vK{uJ) 



e + uj — v + id 
where £(w) = (1 / ir)lmxR(Lu) . Introducing 
1 



€ + uj — v + iS 



and 



7(<) 
7i (t) 



we obtain 

E(c) 



-it/ 2 



(A-2) 

it e^+^+i- 5 )* (A-3) 

du e lult i{uj)n{uj), (A-4) 

1 

duj e iult t(uj)[n(uj) + 1], (A-5) 

dt e ut [(a(t) + (3(t)h(t) 

+a(t)(7i(*)-7(*))]. (A-6) 



But since we can easily prove 7i(t) = 'y(-t), we finally 
obtain eq.(Ell). 
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